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Abstract We have diagnosed the lifetimes of long-lived source gases emitted at the surface and removed in 
the stratosphere using six three-dimensional chemistry-climate models and a two-dimensional model. The 
models all used the same standard photochemical data. We investigate the effect of different definitions of 
lifetimes, including running the models with both mixing ratio (MBC) and flux (FBC) boundary conditions. 
Within the same model, the lifetimes diagnosed by different methods agree very well. Using FBCs versus MBCs 
leads to a different tracer burden as the implied lifetime contained in the MBC value does not necessarily match 
a model's own calculated lifetime. In general, there are much larger differences in the lifetimes calculated by 
different models, the main causes of which are variations in the modeled rates of ascent and horizontal mixing 
in the tropical midlower stratosphere. The model runs have been used to compute instantaneous and steady 
state lifetimes. For chlorofluorocarbons (CFCs) their atmospheric distribution was far from steady state in their 
growth phase through to the 1980s, and the diagnosed instantaneous lifetime is accordingly much longer. 
Following the cessation of emissions, the resulting decay of CFCs is much closer to steady state. For 2100 
conditions the model circulation speeds generally increase, but a thicker ozone layer due to recovery and 
climate change reduces photolysis rates. These effects compensate so the net impact on modeled lifetimes 
is small. For future assessments of stratospheric ozone, use of FBCs would allow a consistent balance between 
rate of CFC removal and model circulation rate. 


1. Introduction 

The abundance of a species in the atmosphere depends on its rate of emission and its rate of chemical (and/or 
physical) loss, i.e., its lifetime. Estimated lifetimes are used to predict future abundances of ozone-depleting 
substances (ODSs) and greenhouse gases (GHGs) and are therefore used to create scenarios used for modeling 
the future atmosphere [e.g., World Meteorological Organization ( WMO ), 201 1 ]. These scenarios provide boundary 
conditions which largely control the time evolution of the atmospheric burden of the source gases, therefore 
constraining projections for the recovery of the ozone layer, for example. Lifetime estimates are also used for 
projections of radiative forcing and climate change. Over the past few years significant doubts have been raised 
over some of the tabulated lifetimes of major ODSs provided in, e.g., World Meteorological Organization/United 
Nations Environment Programme Assessments [e.g., Douglass et al., 2008; Liang et al., 2008]. 

Estimating lifetimes generally requires an atmospheric model. For certain species, with a good measurement 
database, observations can be used to calculate relative lifetimes [e.g., Volk et al., 1997] or to help constrain the 
model. For other species lifetime estimates must be purely model based. In the past two-dimensional (2-D) 
models have been used to estimate the lifetimes of ODSs [e.g., Kaye et al., 1 994; Naik et al., 2000]. Flowever, our 
ability to model atmospheric chemistry and transport has increased greatly over the past couple of decades, 
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Table 1 . List of CCMs and 2-D Model Used in This Study 


Model 


Institutes 


Researchers 


Reference 


and it is important to review and 
update lifetime estimates using 
state-of-the-art models. In particu- 
lar, multidimensional models now 
perform much better in realistically 
simulating the slow stratospheric 
Brewer-Dobson circulation [e.g., 
Stratosphere-troposphere Processes 
and their Role in Climate Chemistry- 
Climate Model Validation Activity 
(SPARC CCMVai), 2010], compared 
to the unrealistically fast circulation 
models used by Kaye et al. [1994], 

It is very important to have the 
best possible estimates of ODS 
lifetimes because of the crucial 
role lifetimes play in guiding 
policymaking in assessing these 
substances. The World Climate 
Research Programme's (WCRP's) 
Stratosphere-troposphere Processes 
and their Role in Climate (SPARC) 
Lifetime Assessment [SPARC, 2013] 
recently conducted a reevaluation 
of the lifetimes of 27 important 
ODSs, their replacements, and GHG 

species based on theory, laboratory studies, observations, and modeling. The species considered were divided 
into those that undergo mainly stratospheric removal (SR) via photolysis or reaction with excited atomic oxygen 
(O^D)) and those that undergo mainly tropospheric removal (TR) via reaction with OH. The 1 1 SR species 
considered comprised chlorofluorocarbons (CFCs), halons, and nitrous oxide, while the 16 TR species included 
hydrochlorofluorocarbons (HCFCs), hydrofluorocarbons (HFCs), and methane. Although some halons can be 
destroyed in the troposphere through photolysis at long-ultraviolet (UV) wavelengths, SPARC [2013] treated 
these as SR tracers due to the nature of their main loss process. SPARC [2013] not only provided a best estimate 
of the lifetimes but also provided a robust estimate of uncertainty, based on the various laboratory, observation, 
and model studies. 


GEOSCCM 

NASA GSFC, USA 

Q. Liang 
A. R. Douglass 
S. Frith 

Pawson et al. [2008] 

GSFC2D 

NASA GSFC, USA 

E. Fleming 
C. Jackman 

Fleming et al. [2011] 

LMDZrepro 3 

CNRS, France 

S. Bekki 
M. Marchand 
J. Burgalat 
D. Cugnet 

Jourdain et al. [2008] 
Marchand et al. [201 2] 

SOCOL 

PMOD-WRC, ETH, 
Switzerland 

E. Rozanov 
A. Stenke 

F. Tummon 

Stenke et al. [2013] 

UMUKCA 

U. Cambridge, UK 
Met Office, UK 
U. Leeds, UK 

P. Braesicke 
N. L. Abraham 
J. A. Pyle 
S. Hardiman 
N. Butchart 
S. Dhomse 

Morgenstern et al. [2009] 
Telford et al. [201 3] 

ULAQ 

U. L'Aquila, Italy 

G. Pitari 

D. lachetti 
G. Di Genova 

E. Mancini 

Pitari et al. [2002] 

WACCM 

NCAR, USA 

D. Kinnison 

Garcia et al. [2007] 


'The LMDZrepro version used here has only 39 vertical levels, the same as 
in the CMIP5 LMDZ climate model [Dufresne et al., 201 3], instead of 50 levels 
as in SPARC [2010]. 


In this paper we discuss the model runs performed for SPARC [201 3]. We present lifetime estimates of 1 1 long- 
lived SR ODSs from six three-dimensional (3-D) Chemistry-Climate Models (CCMs) and one 2-D model from 
groups worldwide for present-day and 2100 conditions (Table 1). The CCMs (or earlier versions thereof) par- 
ticipated in the Chemistry-Climate Model Validation Activity (CCMVai) for the WCRP's SPARC Project [Eyring et al., 
2005], and so the key processes which impact lifetimes (e.g., photolysis and stratospheric circulation) have 
already been evaluated extensively [e.g., Waugh and Eyring, 2008; SPARC CCMVai, 201 0; Morgenstern et al., 201 0; 
Bute hart et al., 2011; Strahan et al., 2011]. 

CCMs participating in recent stratospheric ozone assessment efforts [e.g., WMO, 2011] used prescribed time- 
dependent mixing ratio boundary conditions (MBCs) at the surface to simulate the evolution of ODSs and GHGs 
in the atmosphere. This approach can lead to inconsistencies between the burden of the trace gas and CCM 
transport and loss if these rates differ from the real atmosphere for studies using past observations or to those in 
the model which created a future surface trace gas mixing ratio scenario. We investigate that consistency here 
by using CCMs to investigate differences that arise from using emissions, formulated as flux boundary condi- 
tions (FBCs), to specify the input of ODSs and GHGs. 

The aim of this paper is to complement the updated lifetime estimates for SR species presented in SPARC 
[2013], which are partly based on CCM simulations, by providing further analysis of the models. In particular, 
we investigate the causes of model-model differences and variations in lifetime estimates within the same 
model. We also investigate how lifetimes may change in a future atmosphere. Section 2 describes how the 
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lifetime simulations were formulated. The main results of model evaluation and lifetime estimation are 
contained in section 3. Our results are summarized in section 4. This paper does not attempt to suggest 
recommended lifetimes. For that the reader is referred to SPARC [201 3] which provides best estimates based 
on results of our model runs combined with independent values from observations. 

2. Model Experiments 

We have used results from seven atmospheric models which performed simulations for the recent past, 
present, and future atmosphere. The model experiments are described in this section, along with our 
methods for deriving lifetimes. 

2.1. Participating Models 

We use simulations from six 3-DCCMs and one 2-D model, Goddard Space Flight Center (GSFC2D) (Table 1). All of 
these simulations used the same kinetic and photochemical recommendations from Jet Propulsion Laboratory 
(JPL) 10-6 [Sander eta!., 201 1] and were driven with the same surface boundary conditions for GFIGs and ODSs, 
sea ice, and sea surface temperatures (including GSFC2D which was driven by the zonal mean surface tempera- 
tures from Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM)). Flowever, the models differ in 
dynamical schemes, e.g., advection, convection, cloud parameterization, and the inclusion or treatment of various 
atmospheric processes, e.g., solar cycle, quasi-biennial oscillation (QBO).The six 3-D CCMs all participated in SPARC 
CCMVal-2 and the report from that activity, SPARC CCMVal [2010], contains descriptions of the models and a 
detailed process-based evaluation. Updates to the models for the versions used here are described in SPARC [201 3]. 
In this paper we focus on variations in lifetime estimates due to differences in model dynamics and transport. 

2.2. Lifetime Experiments 

We use a number of model simulations and types of tracers to investigate the relationships between surface 
fluxes, atmospheric burden, removal rates, and lifetime. 

2.2.1. Model Simulations 

We conducted and analyzed the following integrations: One transient simulation from 1960 to 2010 (TRANS) 
and two time slice simulations with 2000 conditions (TS2000) and 21 00 conditions (TS21 00) to calculate present- 
day and future lifetimes, respectively. Simulations TRANS andTS2000 allow us to compare lifetimes calculated in 
a full transient experiment with a steady state experiment. All models used the same mixing ratio boundary 
conditions for long-lived GFIGs from the Coupled Model Intercomparison Project (CMIP) Representative 
Concentration Pathways (RCP) Scenario 4.5 ("historic" scenario before 2000), ODSs according to WMO [201 1], 
and FlFCs based primarily on Velders et at. [2009]. Sea surface temperatures and sea ice concentrations in all six 
3-D CCMs are prescribed as monthly mean boundary conditions following the global sea ice concentration and 
sea surface temperature (FHadlSSTI ) data set provided by the UK Met Office Fladley Centre [Rayner et at., 2003]. 
We present here a brief description of the model simulations; for more detailed information, see SPARC [201 3]. 

TRANS is a 50 year transient run from 1 960 to 2010, based on the definition of the REF-B1 simulation used in 
CCMVal-2 [SPARC CCMVal, 2010, chapter 2; Morgenstem et a!., 2010]. Note that while REF-B1 was a transient 
simulation with mostly stratosphere-only chemistry schemes, the TRANS simulations from four (SOCOL, 
ULAQ, WACCM (Whole Atmosphere Community Climate Model), and UMUKCA) of the six 3-D CCMs were run 
with stratosphere-troposphere coupled chemistry. All forcings in this simulation are taken from observations 
and are mostly identical to those described by Eyring etal. [2006] and Morgenstem etal. [2010] for REF-B1 . This 
transient simulation includes all anthropogenic and natural forcings based on changes in trace gases, solar 
variability, volcanic eruptions, quasi-biennial oscillation (QBO), and sea surface temperatures/sea ice con- 
centrations (SSTs/SICs). Models that do not include a detailed tropospheric chemistry scheme (GEOSCCM, 
GSFC2D, LMDZrepro) prescribed their tropospheric OPI values from the 3-D monthly fields of Spivakovsky 
et at. [2000]. For models that have coupled stratosphere-troposphere chemistry schemes, emissions of ozone 
and aerosol precursors are from the RCP 4.5 scenario [Lamarque et al., 201 1]. 

TS2000 is a 30 year time slice simulation for 2000 conditions, designed to diagnose steady state lifetimes and 
to facilitate the comparison of model output against constituent observations from various measurement 
data sets [SPARC, 2013]. This simulation is conducted with prescribed GFIG, ODS, and FHFC surface boundary 
conditions for 2000, but individual models were run with either repeating or interannually varying solar 
variability, volcanic eruptions, QBO, and SSTs/SICs for 2000 conditions. The final 20 years are used for analysis. 
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TS2100 is similar to TS2000 but for assumed 2100 conditions and used to diagnose steady state lifetimes in a 
future climate with a recovered stratospheric ozone layer and a possibly different Brewer-Dobson circulation. 
This simulation was conducted with prescribed GHG (RCP 4.5 scenario) and ODS surface boundary conditions 
for 2100, but individual models were run with either repeating or interannually varying solar variability, QBO, 
and SSTs/SICs for 2100 conditions. 

2.2.2. Additional ODS Tracers 

Two additional sets of ODS tracers were embedded in the CCM simulations but uncoupled from the full 
chemistry scheme. One set of tracers used realistic surface emission FBCs (FBC), while the other used time- 
independent MBCs (CONST). Although the FBC and CONST tracers are driven with different boundary 
conditions, they are destroyed in the atmosphere with the same kinetics as the corresponding full chemistry 
MBC tracers. Note that potential differences in simulated FBC concentrations do not affect the modeled 
chemistry and physics as these tracers are uncoupled from the interactive chemistry scheme which is based 
on the MBC tracers. 

2 .2 .2.1. FBC Tracers 

We include two stratospheric removal FBC tracers in the CCM simulations for CFC-1 1 (CFC-1 1_FBC) and CFC-12 
(CFC-1 2_FBC), the two most abundant chlorinated ODSs that have ready-to-use bottom-up emission estimates. 
These two FBC tracers are initialized with similar conditions to the MBC tracers at the start of the simulation 
and evolve with geographically resolved surface emission fluxes and atmospheric losses via photolysis and 
reaction with O('D). The model runs also included two tropospheric removal (TR) FBC tracers for CH 3 CCI 3 and 
FICFC-22 which are discussed in SPARC [2013]. 

2 .2 .2.2. CONST Tracers 

Three constant tracers are embedded in the simulations: CFC-1 1_CONST, CFC-1 2_CONST, and N 2 0_CONST. 
These tracers are lost with the same kinetics as the MBC tracers but have a prescribed constant 100 parts per 
trillion by volume (pptv) surface boundary condition. 

We use the combination of MBC, FBC, and CONST tracers from the same model runs to examine how surface 
emission trends, atmospheric abundance, and trace gas distribution impact atmospheric lifetime. 

2.3. Methods for Estimating Lifetimes 

The atmospheric lifetime (r) is generally defined as the ratio of atmospheric burden to the removal rate of the 
species. Therefore, the lifetime of a species depends not only on the properties of the molecule (e.g., ab- 
sorption cross section) but also on its emission history and the atmospheric state. There is no unique value for 
a lifetime of a species and various definitions exist. An instantaneous lifetime is defined at a fixed moment in 
time (or over a short time window) and is estimated from the atmospheric state at that time. In general 
lifetimes based on observations are instantaneous ones valid at the time of the observations. A steady state 
lifetime (r ss ) is the lifetime when the species is at steady state in the atmosphere (i.e., sources = sinks). This 
rarely occurs in the real atmosphere for long-lived source gases, but steady state conditions can easily be set 
up in a model simulation (e.g., in a time slice run). Note that measurements of instantaneous lifetimes can be 
corrected in order to report observational estimates of z ss . Partial lifetimes can also be defined. These can be 
the lifetime of a species against loss by a particular process (e.g., z- phot for loss by photolysis, r 0 iD for loss by 
reaction with O^D), z Q H for loss by reaction with OFI), or the loss of a species within a subregion of the 
atmosphere (e.g., r strat for loss within the stratosphere and r trop for loss within the troposphere). The overall 
lifetime can be calculated from its components by 



^total 


The atmospheric partial lifetime, r atrnos , of a trace gas is calculated in a model using the globally integrated 
sum of its mass-weighted atmospheric loss rate at all locations. Note that except for species that have surface 
losses, e.g., CCI 4 among those considered here, the atmospheric partial lifetime r atmos is equivalent to the 
global atmospheric lifetime z. In this paper all of our quoted lifetimes are atmospheric lifetimes. 

Aside from uncertainty associated with kinetic rates [see SPARC, 2013], uncertainty in lifetime estimates can 
also arise from variations in the representation of model transport and chemical processes. For species that 
are removed mainly by UV photolysis with maximum loss in the stratosphere, lifetimes depend not only on 
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H-121 1 loss rate CFCI 3 loss rate 



0.0 0.8 1.6 2.4 3.2 4.1 0.0 0.4 0.8 1.2 1.6 2.0 

x10‘ 5 molec./cm3/s xIO 3 molec./cm3/s 


CF 2 CI 2 loss rate N 2 0 loss rate 



0.0 0.3 0.6 0.8 1.1 1.4 0.0 1.1 2.2 3.3 4.4 5.5 

xIO' 3 molec./cm3/s xIO’ 1 molec./cm3/s 


Figure 1 . Latitude-pressure cross sections of zonally integrated annual loss rates of Halon-1 21 1 , CFCI 3/ CF 2 CI 2 , and N 2 0 between 2000 and 2005 from the WACCM 
TRANS simulation with warm colors indicating faster loss rates. The solid contours outline the regions within which 95%, 75%, and 50% of the loss occurs. 


the photolysis rate but are also affected by how fast the atmospheric circulation moves air through the 
maximum loss region. 

Because lifetime depends on burden, a potentially important and easily avoided source of model error 
(and model-model differences) relates to the treatment of surface orography. For SR tracers the model burden is 
dominated by the troposphere where the tracer is inert and well mixed. In the extreme case, if a model ignored 
surface orography and assumed a uniform surface pressure of 1013 hPa, it would have an atmospheric mass 
below 400 hPa about 5% larger than a model with realistic topography and a global mean surface pressure of 
around 985 hPa. For such models (here GSFC2D and ULAQ) we can account for this by assuming a constant 
atmospheric mixing ratio from 400 hPa to the surface at 1013 hPa and applying a mass correction factor of 0.95 
to the calculated burden below 400 hPa. For the FBC tracers this issue would manifest itself in a different way. 
Models will have the same emitted tracer mass but in a model which ignores surface orography this mass will be 
distributed over a deeper atmosphere, likely leading to smaller tropospheric mixing ratios. This could lead to a 
smaller concentration flux through the tropopause into the stratospheric sink region and therefore a longer 
lifetime estimate. Flowever, this small potential effect is not easy to treat using archived data and so is ignored 
here. This results in a slight difference in calculated lifetimes between MBC tracers and their corresponding 
FBC tracers. 

3. Results 

To illustrate the importance of chemical loss in different regions, Figure 1 shows the annually averaged 
zonal mean loss rates of four ODSs. CFC-1 1, CFC-12, and N 2 0 are SR species with chemical loss occurring 
solely in the stratosphere. The shortest lived of these, (CFC-1 1), is destroyed lower in the stratosphere 
compared to longer-lived CFC-12 and N 2 0. Although CBrCIF 2 (Halon-1 211) is also primarily destroyed via 
photolysis, the majority of its loss occurs in the troposphere as it can be removed by photolysis at UV 
wavelengths up to 320 nm. For all species considered here, about 75-95% of the atmospheric removal 
occurs between 60°N and 60°S with >50% of the loss occurring in the tropics (within 30° of the equator). 
Therefore, the performance of the models in the tropical region is a key to calculating accurate lifetimes [see 
SPARC, 2013], In particular, the ability of the models to simulate realistic stratospheric tracer transport in the 
tropics is critical and this is investigated in detail in section 3.1. Section 3.2 discusses issues related to the 
quantification of lifetimes from the model runs for present-day conditions. Section 3.3 discusses how lifetimes 
may change in a future climate. 

3.1. Model Mean Age and Comparison With Observations 

Age-of-air is a simple but very useful diagnostic for testing the stratospheric circulation in CCMs and is 
discussed in detail in SPARC CCMVal [2010]. Figures 2 and 3 show the stratospheric modeled age-of-air from 
the CCM runs for present-day conditions and comparison with values derived from balloon and aircraft 
observations [Andrews et ai, 2001 ; Engel et at., 2009]. For the analysis here the modeled age was set to zero at 
the tropical tropopause because the tracers used to diagnose age in these runs were initialized in the lower 
troposphere and therefore include the transit time to the tropopause. This procedure had the largest impact 
on UMUKCA which otherwise produced a tropopause age with a value of around 0.5 years; the shift of 
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GSFC2D ULAQ SOCOL 



GEOSCCM 


WACOM 


UMUKCA 



LMDZrepro 



-50 0 50 


Latitude 


Age of air at 50 hPa 



WACCM 

ULAQ 

GSFC2D 

GEOS 

SOCOL 

LMDZrepro 

UMUKCA 

Obs. 


Figure 2. Mean age-of-air (years) calculated from the last 15 years of theTS2000 runs. Contour interval 0.5 year, (bottom, 
right) The modeled mean age at 50 hPa with derived age-of-air from ER-2 aircraft measurements [Andrews et al., 2001] 
along with their 2 o uncertainty shown as grey shading. 


reference for the other models was smaller than this. This may be partly due to UMUKCA using an altitude- 
based vertical coordinate compared to the rest of the models which use pressure-based coordinates. The 
diagnostics were performed on common monthly mean pressure fields which would have resulted in larger 
interpolation of the UMUKCA fields. In the lower stratosphere (50 hPa) the CCMs tend to agree well with each 
other and with the observed variation in age from the tropics to the poles. One exception is SOCOL, which 
produces mean age-of-air about 1 year younger than the observations at high latitudes. Sensitivity tests with 
a version of SOCOL with higher vertical resolution (90 levels instead of 39) showed that the stratospheric 
circulation is slower and the age-of-air latitudinal gradient increases (not shown), indicating that the young 
age of air in Figure 2 is due at least in part to the coarse-model vertical resolution. LMDZrepro also stands out 
as having air in the lower stratosphere (LS) much older than the other models, although the agreement is 
better higher up and the model produces relatively short lifetimes (see further discussion in section 3.2.4). 
The LMDZrepro version used here has only 39 vertical levels, the same levels as in the CMIP5 Laboratoire de 
Meteorologie Dynamique (LMDZ) climate model [Dufresne et al., 201 3], whereas the version used in SPARC 
[201 0] had 50 levels. The 50-level version did not show this strong bias in the LS mean age-of-air [SPARC, 
2010]. Clearly, the reduction in the vertical resolution in the stratosphere has a large negative impact on the 
quality of the simulated stratospheric circulation. Overall, the models show a large variation in age in the 
middle and upper stratosphere. At around 50 km SOCOL, WACCM, LMDZrepro, GEOSCCM, and UMUKCA 
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Figure 3. Mean age-of-air from present-day CCM simulations for (a) tropics (10°N-10°S) and (b) northern midlatitudes (35°N-45°N). 
Note the different altitude ranges for Figure 3b. The ages are 1 5 year averages of the model zonal mean output from TS2000 
runs. Figures 3a and 3b also show estimates derived from observations [Andrews et al., 2001; Engel etal., 2009] as used in Strahan 
etal. [2011], with their la uncertainty shown as grey shading, (c) The mean age gradient between the northern midlatitudes 
and tropics from the model runs versus observations with their la uncertainties [see SPARC CCMVal, 2010, Figure 5.5]. 


produce ages that are 1 -2 years younger than ULAQ and GSFC2D. The balloon observations used in 
Figure 3 only provide data up to ~30 km but tend to show the models with the older ages are more realistic. 

Figure 3 also shows the gradient in mean age between the northern midlatitudes and the tropics. This diag- 
nostic evaluates the tropical ascent rate [Neu and Plumb, 1999; CCMVal SPARC, 2010]. GSFC2D, ULAQ, and 
UMUKCA show the best agreement with the observations in the 18-27 km (-70-20 hPa) region where CFC-1 1 
experiences significant loss. GEOSCCM and WACCM have gradients that indicate their tropical ascent is some- 
what faster but still within the uncertainty range of the observations over most of this range. SOCOL is the only 
model showing significantly fast ascent throughout the tropical lower and middle stratosphere. LMDZrepro has 
mixed results, with very slow ascent below 1 9 km but fast ascent above. Overall, these results show that ULAQ, 
GEOSCCM, GSFC2D, UMUKCA, and WACCM have reasonable tropical ascent rates. Florizontal mixing (i.e., 
recirculation) into the tropics, which impacts tropical mean age and the ODS burden in the maximum loss region, 
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is not assessed by this diagnostic. SPARC [201 3] has 
further comparison of the model runs with obser- 
vations of some of the target SR species. 

3.2. Model Lifetime Calculations 
3.2.1. Evolution of Lifetimes From 1960s 
to Present 

The time evolutions of the instantaneous lifetimes 
(^transient) of all SR species and halons from the 
TRANS run are shown in Figure 4. Note that while 
the TRANS run extends from 1 960 to 201 0, only 
the results between 1960 and 2006 are shown 
as that is the period for which all model results 
were available. The r transient displays significant 
interannual variations with year-to-year la vari- 
ance of ±3-5%. Overall, the model calculations 
clearly show a decrease in lifetime for all SR spe- 
cies due to the diminishing imbalance between 
surface fluxes and atmospheric losses as the 
atmosphere approaches steady state conditions 
[, Martinerie et ai, 2009]. The contrast in lifetime 
between CFC MBC and CONST tracers, shown in 
Figure 5, clearly demonstrates that the main cause 
of the decrease in CFC-1 1 and CFC-12 lifetimes 
can be attributed to trends in their atmospheric 
concentrations. Secondary causes of the decrease 
in lifetime of the SR species are likely due to the 
combination of higher-altitude 0 3 depletion that 
increases photolytic destruction and any changes 
in the Brewer-Dobson circulation [Fleming et ai, 
2011], The CONST tracers implemented in the 
GEOSCCM and WACCM models suggest that the 
changes in photolysis and circulation together 
explain a decrease of ~5 years (-8% with respect 
to the -57-58 years in 2000s) in CFC-1 1 lifetime 
and -7 years (-7% with respect to the -93-96 years 
in 2000s) in CFC-12 lifetime between the 1960s and 
2000s. Changes in atmospheric concentrations lead 
to a decrease of -15 years in CFC-1 1 lifetime and 
-30 years in CFC-1 2 lifetime over the same period. 
The small difference between the lifetimes of N 2 0 
and its CONST tracer implies that the lifetime of 
N 2 0 is not affected by the relative small change in 
its atmospheric concentration. The -7 year decrease 
in N 2 0 lifetime from 1 960s to 2000s (6% with respect 
to a 115 year lifetime in 2000s) is mainly due to 
changes in photolysis and atmospheric circulation. 


Figure 4. Time evolution of modeled instantaneous 
atmospheric lifetimes (^transient. years) of SR species 
between 1960 and 2006 from the TRANS simulations. 
Model mean lifetimes (thick black lines) and la variance 
(grey shadings) are also shown. The grey-hatched area 
indicates where atmospheric concentrations of the 
corresponding ODSs are close to zero. 
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Figure 5. Time evolution of modeled atmospheric lifetimes of MBC CFC-1 1, CFC-12, and N 2 0 (black lines) between 1960 
and 2010 from the TRANS simulations and their corresponding FBC (blue lines) and CONST (red lines) tracers when available. 
For clarity, the lifetimes are smoothed with a 5 year running mean filter. 


Figure 6 compares model hemispheric mean FBC tracers with observations. The MBC tracers are forced to 
have a tropospheric mixing ratio which matches the observations, while for the FBC tracers this is a predicted 
quantity. Given realistic emission data, the comparison of the FBC tracers with observations, therefore, tests 
the models' ability to reproduce observed loadings through realistic lifetimes and shows the degree of 
model-model variability which would occur if the simulations used FBCs. The FBC tracers display significant 
interhemispheric and seasonal variations. The lifetimes for the FBC tracers are similar to those for the MBC 
tracers (see Table 4 below) because the sink will scale as the atmospheric burden changes so that the lifetime 
just depends on the model transport and chemistry. As the models generally perform well in reproducing the 
observed surface mixing ratios, this shows that their transport and chemical loss rates are reasonable. 
Flowever, despite its relatively long FBC lifetimes, ULAQ simulates a low burden of CFC-1 1 and CFC-12 with 
time which carries on increasing when other models have turned over. Further comparison of the 2000 
steady state and instantaneous lifetimes of MBC, FBC, and CONST tracers is shown in SPARC [201 3], Note that 
the instantaneous lifetimes in 2000 for the listed ODSs are not statistically different from their year 2000 
steady state lifetimes, despite their different trends in atmospheric concentrations (decreasing trends for 
CFC-1 1 , CFC-1 2, CCI 4 , and increasing trends for N 2 0). 

3.2.2. Present-Day Steady State Lifetimes 

Table 2 lists the steady state ODS lifetime for 2000 conditions from the seven models. The models agree fairly 
well in lifetime estimates for many SR species, with calculated lifetimes in general within 10% of the 
multimodel mean. Lifetimes of several major SR species, e.g., CCI 4 , CFC-1 1, CFC-12, N 2 0, from the SOCOL and 
LMDZrepro models are shorter than the other models, most likely due to the fast tropical ascent in the lower 
stratosphere (see sections 3.1 and 3.2.4). 

Table 3 summarizes the present-day multimodel mean steady state r ss and 1 a variance of the 1 1 ODSs, 
along with the final recommendations from SPARC [2013]. The model mean lifetimes and their variances 
(note that these are not uncertainty estimates) are calculated from the inverse of the mean and the 
variance of all available modeled loss rates. Table 3 also lists partial lifetimes in the stratosphere and 
troposphere and partial lifetimes associated with different loss processes, i.e., photolysis and reaction with 
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O('D). Table 4 compares the lifetimes 
diagnosed from different types of model 
tracer. The model mean lifetime esti- 
mates for CFC-11, CCI 4 , and all four 
halons are significantly different from 
the WMO [2011] values. Photolytic 
destruction is the dominant loss process 
for all SR species with reaction with (VD) 
as a minor loss channel for CFCs and 
N 2 0. Reaction with 0( n D) is particularly 
important for CFC-114 and CFC-11 5, 
accounting for 27% and 38% of their 
global loss, respectively. Although 
photolysis is the predominant removal 
process for Flalons, for three (FI-1202, 
FI-121 1, >and FI-2402) of the four studied, 
the majority of the removal occurs in 
the troposphere. (Partial lifetimes with 
respect to different chemical loss pro- 
cesses for all SR species are detailed in 
SPARC [2013, chapter 3]). We note that 
updates to kinetic and photochemical 
parameters in SPARC [2013], which are 
not used here, resulted in significant 
changes in the lifetimes of some of the 
SR species, especially CFC-11 5, Flalon- 
1202, FHalon-1 211, and Flalon-2402. 
These large differences between the 
SPARC final recommendations and our 
model mean lifetimes for these species 
are mainly due to the differences between the recommended reaction rates for CFC-1 15 + O('D) and, 
for the halons, the UV absorption cross sections. The SPARC recommended rates are based on new 
laboratory measurements documented in Baasandorj et al. [2013] and Papanastasiou et al. [2013]. 

3.2.3. Relative Lifetimes From Tracer-Tracer Correlations 

For SR species under steady state conditions tracer-tracer correlations can be used to derive relative lifetimes. 
As explained by Plumb and Ko [1992], two tracers with lifetimes longer than quasi-horizontal mixing time 
scales should be in "slope equilibrium" and produce a compact correlation. Species with lifetimes longer than 
vertical transport time scales will also be in "gradient equilibrium" and the compact correlation will be linear. 



Figure 6. Time evolution of global mean surface concentration (pptv) of 
(top) CFC-1 1 and (bottom) CFC-12 from observations and the modeled 
FBC tracers in the TRANS simulations from 1 960 to 201 0. The observations 
are from stations averaged over the Northern and Southern Flemispheres. 
The model results are also shown averaged over the Northern (solid lines) 
and Southern (dotted lines) Flemispheres. 


Table 2. Present-Day Steady State Atmospheric Lifetimes (r ss , Years) From the TS2000 Runs 


Species 

Chemical Formula 

Model Mean 

GSFC2D 

GEOSCCM 

LMDZrepro 

SOCOL 

ULAQ 

UMUKCA 

WACCM 

CFC-1 1 

CCI 3 F 

55.3 

58.6 

58.3 

49.1 

50.8 

58.6 

56.8 

56.9 

CFC-12 

CCI 2 F 2 

94.7 

103.7 

96.0 

88.2 

84.1 

99.4 

101.5 

93.1 

CFC-1 13 

cci 2 fccif 2 

87.3 

95.4 

88.9 

81.1 

80.5 

92.8 

87.4 

87.4 

CFC-114 

ccif 2 ccif 2 

189 

204 

- 

- 

169 

205 

- 

184 

CFC-1 15 

cf 3 ccif 2 

991 a 

961 

- 

- 

* 

* 

- 

1022 

CCI 4 

CCI 4 

48.6 

50. 7 b 

52.2 

42.0 

41.4 

54.3 

52.0 

51.3 

N 2 0 

n 2 o 

115 

125 

117 

105 

107 

127 

* 

112 

Halon-1202 

CBr 2 F 2 

1.8 

2.1 

- 

- 

- 

- 

- 

1.6 

Halon-1211 

CBrCIF 2 

11.5 

13.5 

- 

- 

- 

9.8 

12.3 

11.0 

Halon-1301 

CBrF 3 

72.2 

77.4 

72.8 

66.0 

67.1 

78.0 

74.7 

71.4 

Halon-2402 

CBrF 2 CBrF 2 

14.5 

13.9 

- 

- 

- 

13.9 

- 

15.9 


a We calculate the mean lifetime of CFC-1 15 as the average of the GSFC2D and WACCM lifetimes as SOCOLand ULAQ specified a too large reaction rate constant 
for O('D) and CFC-1 15 in their calculations. 

‘’The CCI 4 lifetime from GSFC2D is the r phot . There is also a minor loss channel of CCI 4 against OFI with partial lifetime (r OH ) of 161 years [see SPARC, 201 3, chapter 3], 
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Table 3. Present-Day Multimodel Mean Steady State Atmospheric Lifetimes (r ss , Years) Compared to Results From WMO 
[2011] and SPARC [201 3] a 


Species 

WMO [2011] 

SPARC [2013] 

Model Mean 

tss 

Model Variance 
% (Years) 

^strat 

^trop 

^phot 

^OID 

CFC-1 1 

45 

52 

55.3 (7) 

8% (4.2) 

57.0 

1870 

56.4 

2930 

CFC-1 2 

100 

102 

94.7 (7) 

8% (7.3) 

95.5 

11600 

100 

1750 

CFC-1 13 

85 

93 

87.3 (7) 

6% (5.5) 

88.4 

7620 

92.9 

1460 

CFC-1 14 

190 

189 

189 (4) 

10% (18.0) 

191 

19600 

261 

684 

CFC-1 15 

1020 

540 

991 c (2) 

4% (43.1) 

997 

126000 

1590 

2610 

CCI 4 b 

35 

44 

48.6 (7) 

1 2% (5.6) 

50.6 

1230 

48.7 

- 

N 2 0 

1 14 d 

123 

11 5 (6) 

8% (9.0) 

116 

15600 

127 

1180 

Halon-1202 

2.9 

2.5 

1.8 C (2) 

21% (0.4) 

15.3 

2.0 

1.8 

10600 

Halon-1211 

16 

16 

1 1 .5 C (4) 

14% (1.6) 

33.5 

17.3 

11.5 

8040 

Halon-1301 

65 

72 

72.2 (7) 

7% (4.7) 

73.5 

4490 

73.3 

5260 

Halon-2402 

20 

28 

14.5 C (3) 

8% (1.1) 

33.8 

25.1 

14.4 

6790 


The SPARC [2013] lifetimes are partly based on our model runs but also include estimates from observations where 
possible. Also shown are the partial lifetime contributions to the model r ss . Species whose WMO [2011] lifetimes are less 
(or more) than model mean - (+) la variance are highlighted in bold. The number of models which contribute to the 
mean are shown in parentheses in column 4 (see Table 2). 

The CCI 4 lifetime from the models is r atmos , compared with r tota | from WMO [2011] which accounts for ocean and soil losses. 
c Based on JPL-2010 [Sander et al., 2011]. Significantly different lifetimes are obtained using kinetic and photochemical 
parameters updated for SPARC [2013, chapter 3], 
lifetime from Intergovernmental Panel on Climate Change ( IPCQ [2007], 

For SR species the slope of tracer-tracer correlation plot for points below the photochemical loss region 
(i.e., tropical lower stratosphere) can be approximated by the ratio of the tracer lifetimes: 

_ <72 jgi 
T i <7i d<72 

where g, and g 2 are volume mixing ratios. Plumb [1 996] refined this "global mixing" model to account for sub- 
tropical barriers to transport. 

Equation (2) has been applied to derive estimates of relative lifetimes from observations [e.g., Volk et al., 1 997; 
Brown et al., 2013; Laube et al., 201 3]. In these studies the slope of the tracer-tracer correlation is only eval- 
uated in the extratropics to avoid errors related to subtropical transport barriers. Moreover, time trends in the 
tropospheric mixing ratios of the source gases complicate the analysis as the atmosphere will not be in 
steady state. The application of equation (2) to model experiments is more straightforward as we can perform 
idealized steady state experiments with constant source gas surface boundary conditions. 

We have derived relative lifetimes from simulation TS2000 by correlating model tracers with respect to CFC-1 1 
using the last 15 years of the run when the steady state assumption should be valid (see Figure 7). To evaluate 
the lifetime using equation (2) mean tracer values at/below 100 hPa in the extratropics (35°-60° latitude) were 


Table 4. Comparison of Model Steady State Lifetime (r ss , Years) From theTS2000 Simulation With the Transient Lifetime 
(^Transient) for MBC Species and Their FBC and CONST Tracers in Year 2000 From the TRANS Simulation 


Species 


GSFC2D 

GEOSCCM 

LMDZrepro 

SOCOL 

ULAQ 

UMUKCA 

WACCM 

CFC-1 1 

^SS 

58.6 

58.3 

49.1 

50.8 

58.6 

56.8 

56.9 


^transient, MBC 

56.6 

60.0 

47.0 

49.5 

58.2 

56.3 

55.1 


^transient, FBC 

60.0 

58.7 

- 

- 

60.6 

61.4 

55.2 


^transient, CONST 

- 

60.4 

- 

- 

- 

- 

57.1 

CFC-1 2 

r SS 

103.7 

96.0 

88.2 

84.1 

99.4 

101.5 

93.1 


^transient, MBC 

102.1 

98.5 

85.0 

80.5 

101.4 

99.9 

93.3 


^transient, FBC 

99.1 

96.1 

- 

- 

105.0 

109.7 

93.0 


^transient, CONST 

- 

97.5 

- 

- 

- 

- 

92.1 

CCI 4 

r SS 

50.7 

52.2 

42.0 

41.4 

54.3 

52.0 

51.3 


^transient, MBC 

49.2 

53.5 

40.0 

41.0 

53.8 

51.2 

50.1 

N 2 0 

r SS 

125 

117 

105 

107 

127 

- 

112 


^transient, MBC 

124 

120 

- 

105 

130 

- 

113 


^transient, CONST 

- 

119 

- 

- 

- 

- 

112 
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Figure 7. Plots showing the calculation of relative lifetime for CFC-12 from tracer-tracer correlation plots for the seven 
models. Results are shown for average of last 15 years of the TS2000 runs, (columns 1 and 3) The CFC-12 (CF 2 CI 2 ) versus 
CFC-1 1 (CFCI 3 ) correlation plots. The points in orange are fitted to a straight line (equation on top left). For all models except 
UMUKCA the orange points were selected as those between 35°S-60°S and 35°N-60°N, below 100 hPa and within 30% of 
the maximum tracer values in this region. For UMUKCA the orange points were further filtered to use those with the 
minimum value of CFC-12 for any value of CFC-1 1. (columns 2 and 4) The relative lifetime of CFC-12 evaluated using 
equation (2) (see text) and an assumed CFC-1 1 lifetime of 55 years. The lifetime at 1 00 hPa in the tropics (location of blue 
star) is shown in Figure 7 (columns 2 and 4) and used in Table 5. 


used. Table 5 shows the derived relative stratospheric lifetimes, along with the corresponding absolute lifetime 
assuming a r C Fci 3 of 55 years, which is close to the multimodel mean in Table 2. 

The relative lifetimes shown in Table 5 show a generally good level of agreement between the models. For 
example, the relative lifetime of CFC-12 varies from 89.3 years in WACCM to 99.9 years in UMUKCA. This range 
is much narrower than the range of absolute lifetimes of 84.1 to 1 03.7 given in Table 2. This is also the case 
for CFC-1 13 (83.2-88.9 years versus 80.5-95.4 years), CCI 4 (44.6-49.7 years versus 41.4-54.3 years), and N 2 0 
(1 07.5-1 1 7.7 years versus 1 05-1 27 years). For Flalon-1 301 UMUKCA produces a very curved correlation which 
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Table 5. Relative Steady State Stratospheric Lifetimes (r strat ) for Year 2000 Conditions for SR Species From the TS2000 

Model Simulations, Derived From Extratropical Tracer-Tracer Correlations Compared to CFC-1 1 (CFCI 3 ) a 

Species GSFC2D GEOSCCM LMDZrepro SOCOL ULAQ UMUKCA WACCM 


CFC-1 2 
CFC-1 13 
CFC-1 14 
CFC-1 15 
CCI 4 
N 2 0 

Halon-1 211 
Halon-1301 
Halon-2402 


1 .63 (89.9) 
1.51 (83.2) 
3.09 (169.9) 
14.1 (774.4) 
0.88 (48.3) 
1.96 (107.7) 
0.71 (38.8) 

1 .26 (69.2) 
0.67 (36.6) 


1.67(91.9) 

1.55(85.0) 


0.88(48.6) 

2.03(111.8) 

1 .24(68.2) 


1.78(97.9) 
1 .62(88.9) 


0.86(47.6) 

2.14(117.7) 

0.66(36.2) 

1.36(74.7) 


1.64 (90.1) 


0.87 (44.9) 
2.11 (115.8) 

1.31(72.1) 
0.63 (34.5) 


1 .63(89.4) 
1.52(83.8) 
3.20(176.0) 

0.90(49.7) 
2.03 (111.9) 
0.49(27.7) 

1 .29(70.8) 
0.58(31.6) 


1.82(99.9) 
1 .57(86.4) 


0.81(44.6) 

2.10(115.3) 

b 

1.03(56.7) 


1.62(89.3) 

1 .53(84.4) 
3.20(176.2) 

1 8.9(1 042.8) 
0.88(48.5) 
1.95(107.5) 
0.59(32.5) 
1.25(68.5) 
0.65(35.7) 


a Values in parentheses are the corresponding lifetimes assuming r C Fci3 = 55 years. 
b Derived UMUKCA stratospheric lifetime by this method is short (< 1 7 years). 


is hard to fit with an objective straight line. Aside from its low relative lifetime of 56.7 years, for the other 
models the relative lifetimes in Table 5 also show less spread than Table 2 for Halon-1 301 (68.2-74.7 years versus 
66.0-78.0 years). Also, while LMDZrepro, for example, generally produces low absolute lifetimes in Table 2, this 
is not the case for the relative lifetimes. Problems with the representation of ascent and mixing in a model affect 
the lifetime of both species similarly. Thus, the calculation of a relative lifetime reduces the impact of transport 
issues, with the result that relative lifetimes depend more strongly on chemical loss processes. This explains why 
there is better agreement among the model lifetimes in Table 5. Halons 121 1 and 2402 have significant tro- 
pospheric loss; the tracer-tracer correlation method will diagnose stratospheric lifetimes (compare to Table 3). 

The tracer-tracer correlations for UMUKCA in the lower stratosphere, based on the monthly mean output 
were not as compact as the other models and needed to be fitted in a more selective way (see Figure 7 
caption). Note again that UMUKCA is unique among the models presented here in using a hybrid height 
coordinate and a nonhydrostatic dynamical core. This noncompactness in the tracer-tracer correlations may 
also be due to the semi-Lagrangian advection scheme and the associated tracer mass fixer used in UMUKCA. 
The impact of this scheme on advected tracers is discussed by Morgenstern et al. [2009]. 

3.2.4. Variation of Lifetimes With Modeled Age-of-Air 

The lifetime of a SR tracer depends on stratospheric transport, i.e., the Brewer-Dobson circulation. In 
section 3.1 we examined simulated tropical and midlatitude mean profiles and evaluated model tropical 
ascent rate with rates derived from observations. In this section we explore the relationship between the 
derived lifetimes and age-of-air for individual models. 

Figure 8 shows correlation plots of derived lifetimes with the global mass-weighted average of stratospheric 
mean age-of-air (70-1 hPa) for selected SR species. For this analysis the mean age-of-air for each model has 
been reset to zero at the tropical tropopause (see section 3.1 ). Overall, there is a positive correlation: Models 
with older mean age-of-air give longer lifetimes. LMDZrepro tends to be an outlier because its mean age 
below 20 km is too high (Figure 3). For the other models a good straight-line correlation is obtained whether 
the mean age-of-air is averaged from 100 hPa (not shown) or 70 hPa. However, the mass-weighting means 
that the large age-of-air values in LMDZrepro in the lower stratosphere, in a region where SR species are not 
strongly photolyzed, will strongly bias the results. Overall, the good correlation between the modeled life- 
times and mean age-of-air supports the results of relative lifetimes discussed in section 3.2.3 and highlights 
the influence of transport on the modeled lifetimes of SR species. 

Figure 9 shows the correlation between modeled lifetimes and the difference in mean age between two 
levels in the tropical lower stratosphere for CFC-1 1 and CFC-1 2. The chosen altitudes span the main loss 
region for each species. (Although the CFC-1 2 maximum loss region extends to ~5 hPa, we are limited to a 
maximum altitude of 32 km (1 0 hPa) for tropical mean age observations). This quantity is related to the time 
taken for air to ascend through this loss region. Mean age changes in the tropical profile result from vertical 
advection and mixing between the tropics and midlatitudes (i.e., recirculation into the ascending air). The 
stronger correlation seen in Figure 9 than in Figure 8 suggests that lifetimes are particularly sensitive to the 
combined effects of ascent and mixing in the tropical lower midstratosphere. Figures 9a and 9b (solid vertical 
line) represent the observed mean age gradient, with uncertainties shown by the dashed vertical lines. 
Tropical mean age differences from 70 to 20 hPa for the transient experiment for GEOSCCM, UMUKCA, WACCM, 
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GEOSCCM 
— GSFC2D 

— LMDZrepro 
— SOCOL 

— ULAQ 
— UMUKCA 
— WACCM 


Figure 8. Correlation of modeled lifetimes (years) with the global mass-weighted average of stratospheric mean age-of-air 
from 70 hPa to 1 hPa (years) for (a) CCI 4 , (b) CFC-1 13, (c) CFC-1 1, (d) CFC-12, and (e) N 2 0. Mean results are shown from run 
TS2000 (squares) and years 2000-2009 from run TRANS (asterisks). 

and GSFC2D fall within the observational uncertainty. These four models have the most realistic tropical ascent 
over this range (Figure 3) and produce a very small range of CFC-1 1 lifetimes: 56.8-58.6 years from theTS2000 
runs (Figure 9, squares; see Table 2). The most realistic transport over the 70-10 hPa range is seen in GEOSCCM, 
WACCM, and GSFC2D, with ULAQ and UMUKCA just outside the uncertainty range. The CFC-1 2 lifetimes 
calculated from these five TS2000 model runs span 93.1-103.7 years. The models with tropical ascent rates 
that are too fast (SOCOL and LMDZrepro) have the shortest lifetimes for CFC-1 1 (<51 years) and CFC-12 
(<89 years) and have tropical mean age gradients that fall outside of the observed range. Given that the LMDZ 
lifetimes are lower than most other models in spite of relatively old air below 1 9 km (70 hPa), we conclude that 
fast ascent through the key loss region (i.e., above 70 hPa) is the primary factor influencing lifetimes. 

3.3. Future Lifetimes Estimates 

Atmospheric changes could modify tracer lifetimes in various ways in the future atmosphere. For example, a 
speedup in the Brewer-Dobson circulation in the critical region for stratospheric loss would be expected to 
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Figure 9. Correlation of modeled lifetimes (years) with the differ- 
ence in mean age-of-air between two levels in the tropics for (a) 
CFC-1 1 (age difference between 70 hPa and 20 hPa), and (b) CFC-12 
(age difference between 70 hPa and 10 hPa). Results are shown 
from run TS2000 (square) and years 2000-2009 from run TRANS 
(asterisk). Also shown are estimates of the observed tropical mean 
age difference (vertical solid line) and its la uncertainty (dashed 
vertical lines) using the data shown in Figure 3a interpolated to 10, 
20, and 70 hPa. 


lead to a reduction in the lifetime of photo- 
lytically removed species [ Butchart and 
Scaife, 2001]. Lin and Fu [2013] analyzed 
the Brewer-Dobson circulation changes in 
CCMVal-2 model experiments and found 
that while models consistently simulated an 
acceleration of both the shallow and deep 
branches of the Brewer-Dobson circulation 
(70-30 hPa and above 30 hPa), the acceler- 
ation of the deep branch was much smaller. 
They also noted significant differences 
between models. 

Table 6 compares the lifetimes for TS2000 
and TS2100 for six models that conducted 
both runs. The models predict different 
responses in future atmospheric lifetimes 
of SR species. GEOSCCM, SOCOL, and 
UMUKCA show a decrease in lifetimes in 
most of the SR species, while WACCM 
shows a slight increase in all SR species. 
The future lifetime changes predicted by 
GSFC2D and ULAQ are somewhat mixed, 
with increases for some and decreases for 
others. A comparison of change in mean 
age-of-air between the simulations for 
2100 and 2000 suggests that all models, 
except WACCM, show a speeding up of 
the Brewer-Dobson circulation although 
the magnitude of this change varies 
(Figure 1 0). GEOSCCM, GSFC2D, and SOCOL 
show only a small decrease in mean age 
(~1 —3 months), while UMUKCA and ULAQ 
show larger changes. WACCM shows a 
slight speedup (—0—1 month) of the shallow 
branch of the Brewer-Dobson circulation 
but a slowdown (—0—2 months) of the deep 
branch. Flowever, such a weak change in 


Table 6. Steady State Atmospheric Lifetimes (r ss , Years) for 2000 and 21 00 Conditions From the TS2000 and TS21 00 Simulations From Six Models 

Model Mean GSFC2D GEOSCCM SOCOL ULAQ UMUKCA WACCM 


Species 

2000 a 

2100 

2000 

2100 

2000 

2100 

2000 

2100 

2000 

2100 

2000 

2100 

2000 

2100 

CFC-1 1 

56.5 

55.2 

58.6 

57.7 

58.3 

54.7 

50.8 

49.3 

58.6 

57.6 

56.8 

55.2 

56.9 

57.7 

CFC-12 

95.8 

95.3 

103.7 

103.5 

96.0 

90.6 

84.1 

81.9 

99.4 

104.2 

101.5 

98.8 

93.1 

96.7 

CFC-1 13 

88.5 

87.8 

95.4 

94.4 

88.9 

83.4 

80.5 

79.1 

92.8 

96.6 

87.4 

85.5 

87.4 

90.6 

CFC-1 14 

189 

193 

204 

206 

- 

- 

169 

171 

205 

205 

- 

- 

184 

192 

CFC-1 15 

991 

1015 

961 

981 

- 

- 

- 

- 

- 

- 

- 

- 

1022 

1052 

CCI 4 

49.9 

48.6 

50.7 

49.9 

52.2 

48.7 

41.4 

40.2 

54.3 

52.9 

52.0 

50.4 

51.3 

51.8 

N 2 0 

117 

118 

125 

125 

117 

111 

107 

107 

127 

132 

- 

- 

112 

118 

HI 202 

1.8 

2.0 

2.1 

2.5 

- 

- 

- 

- 

- 

- 

- 

- 

1.6 

1.6 

HI 21 1 

11.5 

11.9 

13.5 

13.4 

- 

- 

- 

- 

9.8 

10.0 

12.3 

13.5 

11.0 

11.4 

HI 301 

73.4 

73.1 

77.4 

76.8 

72.8 

68.6 

67.1 

66.1 

78.0 

79.2 

74.7 

76.2 

71.4 

73.5 

H2402 

14.5 

14.8 

13.9 

14.0 

- 

- 

- 

- 

13.9 

14.2 

- 

- 

15.9 

16.4 


a Model mean lifetimes for 2000 in this table are calculated without LMDZrepro results for a fair comparison with the 2100 lifetimes. 
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GSFC2D: dAge (21 00-2000) GEOSCCM : dAge (21 00-2000) 



SOCOL: dAge (2100-2000) 



UMUKCA: dAge (2100-2000) 



Figure 10. Difference in mean age-of-air (months) between 
perform TS2100 run.) 



ULAQ: dAge (2100-2000) 



WACOM: dAge (2100-2000) 



CCM runs for 2100 and 2000. (Note: LMDZrepro did not 


circulation in WACCM might not be statistically significant. On the other hand, all models predict a 
decrease in photolysis (J) rates under 2100 conditions with the recovered ozone layer and a cooler 
stratosphere, though the relative changes in J rates vary between models and species. ULAQ and UMUKCA 
predict the strongest reduction in J rates, e.g., -13% decrease in -/cfc-ii in the tropical lower stratosphere 
between 1 0 and 1 00 hPa with respect to the 2000 conditions, while GSFC2D, WACCM, and GEOSCCM show 
a much weaker reduction of 5-6%. SOCOL displays a medium response of -9% reduction in the corre- 
sponding J CFC . ql rate. From a multimodel mean perspective, the decrease in lifetime due to a faster Brewer- 
Dobson circulation cancels with the increase in lifetime due to weaker photolysis. Therefore, the modeled 
atmospheric lifetimes for SR species under 2100 conditions do not change significantly from the lifetimes 
under 2000 conditions. 

It is important to note that the projections of ODS lifetime changes in a future climate depend upon 
the choice of the GHG scenario used to drive the TS2100 run. We used RCP 4.5 but higher levels 
of CH 4 , N 2 0, and C0 2 from RCP 6.0 (the more likely scenario) or RCP 8.5 will lead to a different 
response in atmospheric circulation and ozone concentration, therefore a different future lifetime for 
the targeted ODSs. 
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4. Conclusions 

We have used results from seven global models (six 3-D CCMs and a 2-D model), integrated with the same 
standard photochemical data to estimate atmospheric lifetimes of source gases emitted at the surface and 
removed in the stratosphere (SR tracers). We diagnosed both instantaneous and steady state lifetimes for 
recent past, present-day, and 2100 conditions. Analysis of these model simulations has identified the key 
processes which control the lifetime estimates of SR tracers in models. 

SPARC [2013] has used these model results, in combination with observations and updates to photo- 
chemical data, to produce recommended lifetimes (and uncertainties) for SR and tropospheric removal 
(TR) tracers. The methods used to derive the best estimates of lifetimes are described in the SPARC 
report. In this paper we have analyzed model experiments to understand the controlling factors for lifetime 
estimates and to understand the behavior of different types of tracer within the models. In both our 
model analysis and the final recommendation in SPARC [2013], the atmospheric lifetimes of the major 
Cl-containing ODSs considered (CFC-1 1, CFC-1 2, CFC-1 1 3, and CCI 4 ) all increased compared to WMO [201 1], 
albeit only marginally for CFC-1 2. This implies that the predicted rates of stratospheric chlorine decay and 
0 3 recovery will decrease. 

The transient lifetimes of, for example, CFC-1 1 and CFC-1 2 before the 1990s are larger than steady state 
lifetimes but decrease to approach the steady state lifetime after this time. N 2 0 shows little differences be- 
tween transient and steady state lifetimes. The strong increase in CFC emissions over the period before the 
1990s causes this difference when the atmosphere was far from steady state. During the decay phase of CFCs, 
their long atmospheric lifetimes mean that the departure from steady state is much less, and so the transient 
lifetime remains a good approximation for CFCs even after emissions start to decrease. 

Experiments which compared lifetimes derived from flux boundary condition (FBC) tracers to those derived 
from mixing ratio boundary condition (MBC) tracers within the same model gave very similar results. This can 
be explained because, at steady state, the lifetime depends strongly on the modeled chemical loss processes 
which are the same in both cases. The consequence of running with FBCs is that the models produce 
different tracer burdens which scale with the model-model difference in lifetimes. 

Use of FBCs allows a consistent balance between rate of source gas removal and model circulation 
rate and would present an improvement for future CCM studies of ozone recovery. With FBCs the 
evolution of chlorine loading, for example, will be consistent with the dynamics and chemistry 
which also affects ozone. Flowever, the implementation of FBCs is not necessarily straightforward. 
Realistic emission estimates are required for all FBC species, along with accurate quantification of any 
tropospheric sinks (e.g., for CCI 4 ). Furthermore, use of FBCs will lead to larger model-model spread in 
predicted ozone recovery due to the different time scales of source gas removal. Although that will 
be a reflection of real model-model uncertainty, it could complicate the identification of other model 
errors or uncertainties. 

Estimated lifetimes of stratospheric removal species show a range of values between models, but there is a 
correlation between the strength of the modeled Brewer-Dobson circulation (as given by mean age-of-air) 
and lifetime. In general, younger mean age correlates with a shorter lifetime, with the tropical mean age 
difference over the altitude range where maximum loss occurs being most critical. 

The model comparisons have revealed a minor but important issue concerning the accurate estimate of 
lifetimes. The modeled burden of an SR tracer is affected by the assumed surface orography. A model which 
assumes a globally uniform surface pressure of 1013 hPa has more mass and therefore more well-mixed SR 
tracer in the lower troposphere than a model with a realistic mean surface pressure of 985 hPa. This feeds 
directly into the lifetime estimate in the calculation of burden/loss rate. 

For simulations with one possible representation of a 2100 atmosphere (RCP 4.5) a thicker (recovered) ozone 
layer in 2100 leads to reduced photolysis which in itself would lead to longer lifetimes. Flowever, most models 
show younger mean age throughout the stratosphere, suggesting a possible speedup of the Brewer-Dobson 
circulation. For other models the circulation change is not so clear. These two processes of reduced photolysis 
and circulation changes lead to canceling effects in the models, and hence, there is no clear trend in the 
lifetimes of stratospheric removal tracers between 2000 and 2100. Flowever, these impacts will depend on 
the future GFIG scenario used. 
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